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Abstract 

We use the functional renormalization group to analyze the temperature dependence of the 
quasi-particle scattering rates in the two-dimensional Hubbard model below half-filling. Using a 
band structure appropriate to overdoped Tl2Ba2Cu06+a; we find a strongly angle dependent term 
linearly dependent on temperature which derives from an increasing scattering vertex as the energy 
scale is lowered. This behavior agrees with recent experiments and confirms earlier conclusions on 
the origin of the breakdown of the Landau Fermi liquid near the onset of superconductivity. 

PACS numbers: 71.10.Hf, 71.27.+a, 74.72.-h 
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Recently Abdel-Jawad and collaborators [l| reported a striking correlation between charge 
transport and superconductivity in heavily overdoped high-temperature superconducting 
cuprates. They found that the onset of superconductivity with doping coincided with the 
appearance in the normal phase of strong anisotropic quasiparticle scattering. The applica- 
tion of a magnetic field to suppress the superconductivity revealed that the anisotropic term 
in the inplane transport scattering rate was linear in temperature violating the perturbative 
quadratic dependence characteristic of a Landau Fermi liquid. Earlier investigations of a 
2-dimensional Hubbard model on a square lattice using a functional renormalization group 
(RG) method found that rf-wave pairing in the overdoped region of the phase diagram was 



driven by the appearance at low energies and temperatures o: 
ing vertex in the particle-particle and particle-hole channels 



'_a. strongly anisotropic scatter- 
3|, IJ] . Further investigations 



revealed that the self-energy is also anisotropic js, [g], Q, [sl. In this letter we present an ex- 
tensive RG study of the doping and temperature dependence of the quasi-particle scattering 
rate with the pairing instability suppressed, and compare our results with the experiments of 
Ref . [l| . We find a rising anisotropic scattering vertex that gives an anisotropic contribution 
to the decay rate with a linear, not quadratic, temperature dependence when the pairing 
divergence is suppressed. We wish to stress that this breakdown of standard Landau Fermi 
liquid behavior is not associated with a divergent density of states from a van Hove singular- 
ity at the Fermi energy but instead is due to strong scattering processes at large momentum 
transfer which appear in the RG flows as a precursor to the Mott insulating behavior at 
half filling. The origin of anomalous decay rate in our calculations is quite different from 
that proposed by Metzner and co-workers [91] who put forward a model based on small angle 
scattering near to a Pomeranchuk instability. 

The experiments by Abdel-Jawad and collaborators jl| were carried on well characterized 
Tl2Ba2Cu06+x samples. Their studies of the interlayer angle-dependent magnetoresistance 
(ADMR) provided also detailed Fermi surface information which we use to fix the band 
parameters in the RG calculations. We use a moderate starting value of the onsite repulsion, 
U, in the one loop RG equations. Our goal is a qualitative rather than a quantitative 
description which would require a larger value of U and multi-loop corrections to the RG 
flow equations. In the experiments a magnetic field is applied to suppress superconductivity 
and study the normal phase down to low temperatures. However it is technically very 
difficult to introduce a magnetic field into the RG calculations. Instead we introduce an 
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elastic scattering term to suppress the (i-wave pairing instability in the calculation of the 
RG flows of the scattering vertex. This vertex is then used as input into a standard lowest 
order calculation of the quasiparticle decay rate. 

The kinetic energy applicable to Tl2Ba2Cu06+x takes the form 

e{kx,ky) = —2ti {cos kx + cos ky) + 4:t2 {cos kx cos ky) 

+2t3 (cos 2kx + cos 2ky) + 4t4(cos 2kx cos ky 

+ cos 2ky cos kx) + 4^5 (cos 2kx cos 2ky) , (1) 

with ti = 0.181, t2 = 0.075, = 0.004, U = -0.010, and = 0.0013(eV). 



Our approach relies on the functional RG equation for the one-particle irreducible (IPI) 

n □ 

generating functional r[$], which is derived in [J, 110|, and which leads to a hierarchy of 
coupled flow equations for the IPl vertices after a suitable expansion of the functional. 
We use a Wilsonian flow scheme with a sharp momentum cutoff. It was shown later that 



this method underestimates effects of small wavelength scattering [ll|], but also that these 
processes are important only if the Fermi surface is close to a van Hove singularity, which 
is not the case in the doping regime studied here. In order to solve the flow equations, the 
hierarchy has to be truncated, and in the following we will use the standard truncation of 
neglecting all vertices with more than four legs. In this approximation, the only quantities 
appearing in the calculation are the self-energy T,\{k) and the 4-point vertices Va(A;i, k2, ks). 
The ki also contain the frequency, ki = {uJi,]<.i). All propagators contain a sharp infrared 
cutoff XA(k) = 6'(|e(k)| — A) in momentum space, where the flow parameter A flows from 
A = cxD to A = with the initial condition Voo{ki, k2, k^) = U . As these equations are still 
too complicated to be solved, we introduce some further approximations by neglecting the 
frequency dependence of all vertices and by discretizing their momentum dependence. The 
latter is done by dividing the Brillouin zone into elongated patches each of which contains 
a part of the Fermi surface. The momentum dependence of the vertices is approximated by 
a step function which is constant in each patch. The vertices are calculated at a reference 
point in each patch, which we choose to lie where the Fermi surface crosses the center of the 
patch, as shown in FIG. [Ufor the case of a hole concentration p = 0.30. 



Typically, at low energies the truncated flows diverge at some finite energy scale A. The 
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FIG. 1: (Color online) Fermi surface (solid line) and discretization of the BZ for p = 0.22. The 
boundaries of the patches (labelled by 1,2,..., 40) are indicated by the dashed lines. All vertices 
are evaluated at the points marked by the red dots, and are taken to be constant within each patch. 



leading divergence can be interpreted as the dominant instabi 



ity, and the scale at which 



In the regime of interest 



the divergence occurs gives an estimate of the corresponding Tc 
here, d-wave superconductivity is the leading instability, and in our approximation takes 
the values = 0.26ti for p = 0.15, = 0.22ti for p = 0.22, and = 0.16ti for p = 0.30. 
These temperatures are way too high, mainly because we neglect self-energy corrections 
in the flow of the scattering vertex. Nevertheless, grows with decreasing hole doping 
reproducing qualitatively the experimental results [l|. 

The experiments were carried out in a high magnetic field, which suppresses superconduc- 
tivity and allows to access the normal state down to low temperatures. However, including 
a magnetic field into our RG calculation is very difficult, so that we choose a different way to 
suppress superconductivity, namely we introduce an isotropic scattering rate I/tq into the 
free part of the action. This smears out the Fermi distribution at the Fermi surface, which 
in turn regularizes the loop integrals and subsequently the flow of the 4-point vertices. This 
scattering rate will only be included in the flow equation for the 4-point vertex, whereas the 
flow equation for the self-energy is left unaltered. We found that for our choice of f/ = Ati 
and the range of temperatures (T > 0.004ti) and dopings {p > 0.15), a scattering rate of 
I/tq = OMi is sufficient to suppress the divergences associated with superconductivity, so 
that on average the vertices are comparable to the bandwidth and the largest vertices do 
not grow larger than ^ 3 x bandwidth (Fig. [2]). 

In FIG.|2]a typical result of our calculations for the 4-point vertex VA(ki, k2, ka) at energy 
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FIG. 2: (Color online) Characteristic momentum dependence of the renormalized vertex 
FA(ki,k2,k3)/ti for p = 0.22, T = OMh, I/tq = 0.2ti at A = 0. In the figure, the depen- 
dence on the two ingoing wave vectors (ki,k2) is shown, where the outgoing wave vector ks is 
taken to lie in patch 1 close to (vr, 0) (cf. Fig. [T|) and k4 is fixed by momentum conservation. 

scale A is displayed for a fixed outgoing wavevector ks close to (vr, 0) as a function of the 
two incoming wavevectors (ki,k2). The remaining outgoing wavevector is determined by 
momentum conservation allowing for umklapp processes. The strongest scattering processes 
occur for a momentum change of (vr, tt). These scattering vertices are approximately oc A~^/^ 
for small values of the cutoff A ~ T. However, for very small cutoff, the square root 
divergence is suppressed due to the presence of the isotropic scattering rate I/tq. 

As we are only interested in the scattering rates at the Fermi surface, which are given 
by Q'S(k G FS,uj — > + i6), we will restrict the calculation of the self-energy to this 
quantity in the following. Obviously, the frequency-dependence of S cannot be neglected in 
the calculation. On the other hand, if we neglect the frequency-dependence of the 4-point 
vertices, it is clear from the structure of the flow equations that S will also be frequency- 
independent, as only Hartree and Fock diagrams are included. However, as shown earlier 
by one of usj^, one can overcome this difficulty by replacing the 4-point vertex appearing 
in the self-energy flow equation by the integrated flow equation of the vertex, schematically. 
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(2) 



where in our approximation the single-scale propagator [4, llOl and the full propagator 
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FIG. 3: Two-loop diagrams contributing to the self-energy. Only diagrams a) and b) contribute to 
the scattering rates. 

G\ are related to the free propagator Go by 

'S'a = XaGo and Ga = Xa^o, (3) 

respectively. The RHS of eq. ([2]) depends on A only through the cutoff xa- After a 
partial integration with respect to A and after explicitly inserting a sharp cutoff XaO^) = 
e(|e(k)| - A) we have 

Sa=o = J rfA^(|e(ki)| - A)5(|e(k2)| - A)^(A - |e(k3)|) 

xVlGo{h)Go{h)Go{h), (4) 

and summation and integration over internal momenta and Matsubara frequencies is implied. 
Thus the first propagator has support above, the second at, and the third below the cutoff 
A. 

The diagrams corresponding to this equation for the self-energy are shown in FIG. [3l 
As we are interested in the scattering rates at the Fermi surface, we need only consider 
diagrams a) and b), because the contribution of diagram c) is real for external frequencies 
uj + i6. For a) and b), for external frequency lim^^^QUJ + iS, we obtain an imaginary part 
oc 6 (e(k3) — e(k2) — e(ki)), reflecting energy conservation. 

It turns out that neglecting the flow of the 4-point vertices, i.e. setting Va = [/ in eq. (El), 
is equivalent to a second order perturbative calculation of the scattering rate, which gives 
a behavior away from van Hove singularities. All deviations from the Landau theory 
scaling form may be attributed to the renormalization of the 4-point vertices. 

Based on eq. (jl]) we calculate both the temperature and the doping dependence of the 
angle-resolved quasi-particle scattering rates at the Fermi surface using the method explained 
above. We find that the scattering rates are anisotropic for all choices of parameters. The 
precise shape of the angle dependence changes with doping, but does not change very much 
with temperature, as shown in FIG. |H In general, we find that in the nodal direction 
(0 = 7r/4) the scattering rates have a minimum, and increase towards the anti-nodal direction 
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FIG. 4: (Color online) Angular dependence of the anisotropic component of the quasi-particle 
scattering rate on segment of the Fermi surface for p = 0.15 (dashed line) and p = 0.30 (solid line). 
The scattering rates are normalized to unity in the anti-nodal direction. The dots are the values 
at different temperatures, the lines are the temperature averaged angular dependence. 

(0 = 0). The size of the anisotropy grows as doping is decreased. Calculations without the 
background scattering rates as a regulator showed that also Tc increases with lower doping. 
This is in accord with the results of the magnetoresistance measurements carried out by 
Abdel-Jawad et al.[l|, where it was found that with decreasing doping both T^. and the 
anisotropic part of the scattering rates increase, whereas the uniform component remains 
constant. 

Separating the scattering rates into an isotropic and an anisotropic part, we write 



where 1/rj = min0l/r($) so that l/ra(0) > 0. We characterize the T-dependence of the 
anisotropic part by averaging over the angle. 



which makes sense as the angular dependence of the anisotropic part is approximately in- 
dependent of temperature (FIG.H]). Using these definitions, we find that the T-dependence 
of I/tj and (I/to) can be fitted very well by a quadratic polynomial, as shown in FIG. |5l 
In the same figure, one sees that for p < 0.25, {I/to) becomes linear in T with a coeffi- 
cient which increases with decreasing hole doping, whereas the isotropic part 1/rj is always 
dominated by a quadratic term which does not change much with doping. Thus our calcu- 
lations reproduce the main features of the striking correlation between charge transport and 
superconductivity reported by Abdel-Jawad and coworkers. 
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FIG. 5: (Color online) Temperature dependence of (a) the anisotropic and (b) the isotropic com- 
ponent of the quasi particle scattering rate at the Fermi surface for different values of the hole 
doping. The solid lines are fits to quadratic polynomials in T. 

Abdel-Jawad et al. [l| also commented on the contrast between the transport properties 
of their overdoped samples which show the shortest lifetimes in the antinodal directions and 
the low temperature ARPES data showing shortest lifetimes in the nodal directions in the 
superconducting phase of overdoped Tl2Ba2Cu06+x samples jl2l|. However, Wakabayashi et 



al. 



131 ] have attributed this latter behavior to a reduction in the elastic scattering rate 



near the antinodes due to a coherence effect which acts strongest at the (superconducting) 
gap energy. Multiple scattering in this case leads to a depression of the elastic scattering 
rate from impurities, defects etc. Such processes are quite distinct to the intrinsic inelastic 
scattering processes in the normal state discussed here. 

These results signal a clear breakdown in Landau-Fermi liquid behavior which leads to a 
universal dependence. We also wish to stress that this result is not due to a proximity to 
the van Hove singularity at the saddlepoint of the band structure. This lies well below the 
Fermi energy at energies ^ T. Further, the increase in the linear term in T with decreasing 
hole dopings occurs as the energy of the van Hove singularity moves further away from the 
Fermi energy. In our calculations the breakdown of Landau theory arises from the increase 
in the 4-point vertex with decreasing energy scale. Note that this increase is not restricted 
to the d-wave pairing Cooper channel since the divergence in this channel is suppressed in 
our calculations. Examination of the RG flows shows that several channels in the 4-point 
vertex grow simultaneously, e.g. particle-hole and particle-particle umklapp processes. This 
phenomenon is not simply a precursor of d-wave superconductivity but rather signals that 
a crossover to strong coupling in several channels of the 4-point vertex is responsible for 
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the breakdown of the Landau-Fermi hquid behavior. This simultaneous enhancement of 
several channels through mutual reinforcement was earlier identified as a key feature of the 
anomalous Fermi liquid in the cuprates and associated with the onset of resonant valence 
bond (RVB) behavior 0, [l^. For the parameters used above, this crossover takes place at 
p ^ 0.14. However, due to the rich momentum dependence of the 4-point vertex, the nature 
of the strong coupling phase is unclear and requires further investigations. 

In conclusion, the RG calculations presented here demonstrate that the anomalous be- 
havior of the inplane quasi-particle scattering rate revealed by the ADMR experiments [l| on 
overdoped cuprates can be found as an intrinsic feature of the doped Hubbard model already 
at weaker interaction strengths. Likewise, the positive correlation of the critical temperature 
for c/-wave superconductivity with the strength of the scattering rate anisotropy comes out 
of the RG treatment without additional assumptions. 

We are grateful to N. Hussey, L. Taillefer, A. Katanin, and W. Hanke for useful discus- 
sions. This study was financially supported by the Swiss National fonds through the NCCR 
MaNEP. CH acknowledges financial support by the DFG research unit 538. 



[1] M. Abdel-Jawad et al., Nature Physics 2, 821 (2006). M. Abdel-Jawad et al., Phys. Rev. Lett. 

99, 107002 (2007). 

[2] D. Zanchi, H. J. Schulz, Europhys. Lett. 44, 235 (1997). 

[3] C. J. Halboth, W. Metzner, Phys. Rev. B 61, 7364 (2000); Phys. Rev. Lett. 85, 5162 (2000). 

[4] C. Honerkamp et al., Phys. Rev. B 63, 035109 (2001). 

[5] D. Zanchi, Europhys. Lett. 55, 376 (2001). 

[6] C. Honerkamp, Eur. Phys. J. B 21, 81 (2001). 

[7] A. A. Katanin, A. P. Kampf, Phys. Rev. Lett. 93, 106406 (2004). 

[8] D. Rohe, W. Metzner, Phys. Rev. B 71, 115116 (2005). 

[9] W. Metzner et al., Phys. Rev. Lett. 91, 066402 (2003); L. Dell'Anna, W. Metzner, Phys. Rev. 
Lett. 98, 136402 (2007). 

[10] C. Honerkamp, M. Salmhofer, Progress in Theoretical Physics 105, 1 (2001). 

[11] C. Honerkamp, M. Salmhofer, Phys. Rev. B 64, 184516 (2001). 

[12] M. Plate et al., Phys. Rev. Lett. 95, 077001 (2005). 



9 



[13] K. Wakabayashi et al., Phys. Rev. B 72, 214517 (2005). 
[14] A. Lauchli et al., Phys. Rev. Lett. 92, 037006 (2004) 



10 



